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^ ■ Abstract 

The paper contains a systematic, model-independent treatment of electric dipole (El) transitions 

in heavy quarkonium. Within the effective field theory framework of potential non-relativistic 

T^ ', QCD (pNRQCD), we derive the complete set of relativistic corrections of relative order v'^ both 

Ph. 
f->i ' for weakly and strongly-coupled quarkonia. The result supports and complements former results 

D ■ 

r~| ' from potential model calculations. 
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I. INTRODUCTION 

In recent years, BES, the B- factories and CLEO have improved on almost all heavy- 
quarkonium radiative transition modes and measured many of them for the first time. Part 
of this impressive progress can be found collected in two reviews by the Quarkonium Working 
Group [ll, I2I]. Among the most recent developments we mention the first measurements 
of K -> Vol by BES [3], T(25) -^ Vbl hy BABAR [4], T(35) -^ Vbl hy CLEO |5|, 
hb ^ Tjb'-f by BABAR and BELLE [6], and the measurements of the branching fractions for 
electromagnetic transitions of the Xb states by CLEO |7| and BABAR |8|. 

Electric dipole (El) transitions are transitions that change the orbital angular momentum 
of the state by one unit, but not its spin. Hence, the final state has a different parity and 
C parity than the initial one. Typical El transitions are n'^^Pj — t- n'^Si 7, e.g. Xcj -^ J/i^l 
or Xbj -^ '^(l'S')7, and n'^Pi — > n^S'07, e.g. he — )■ rjc'j or hh — > rjh'y. El transitions 
happen more frequently than magnetic dipole (Ml) transitions. The branching fraction of 
El transitions can indeed be quite significant for some states, e.g. the branching fraction for 
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Xbi -^ T(15) 7 is (35 ± 8)% and the branching fraction for Xb2 ^ T(1S) 7 is (22 ± 4)% |9|. 
El transitions depend on the wave function already at leading order and therefore provide 
a direct insight into the quarkonium state. 

Electromagnetic transitions have been treated for a long time by means of potential mod- 
els using non-relativistic reductions of phenomenological interactions (see e.g. [10], which 
will be our reference work in the field). A recent comprehensive analysis can be found in jllj. 
However the prog ress made in effective field theories (EFTs) for non-relativistic systems, like 



quarkonium, 12| , and the new large set of accurate data ask for model- independent analy- 



ses. Ml transitions have been treated in an EFT framework in 13|. In this work, we extend 
that previous analysis to El transitions, for which a model-independent treatment has been 
missing so far. 

Effective field theories for quarkonium radiative transitions are built on a hierarchy of 
energy scales [13|: the heavy-quark mass m, the relative momentum of the bound state 
p ~ mv and the binding energy E ~ mv'^, where f ^ 1 is the heavy-quark velocity in the 
center of mass frame. The condition f ^ 1, which works better for bottomonium (f^ ^ 0.1) 
than for charmonium (f^ ^ 0.3), implies m ^ p ^ E. For transitions that involve a 
change in the principal quantum number, the photon energy, k^, scales like k^ ~ mv"^. This 



counting will be assumed throughout the paper, although, for transitions between states with 
the same principal quantum number, the photon energy is smaller. A smaller photon energy 
may be implemented easily in the final expressions of the transition widths by suppressing 
terms proportional to k^ accordingly. Observables like the transition widths are organized 
in an expansion in f ^; the main purpose of the paper is to provide, in an EFT framework, 
an expression for El transition widths valid up to relative order v"^. The work is partially 
based on jl4| . 

The paper is organized as follows. In Sec. |lll we construct the EFT Lagrangian relevant 
for El decays up to relative order f^. The matching is performed in Sec. Illli In Sec. IIVI 
we set up the calculation of the transition rates. In Sec. |Vl we calculate the relativistic 



corrections and, in Sec 



VTI we write the widths up to order f^. There we also compare our 



results with those in 10|. Finally, in Sec. IVIIt we conclude and discuss future applications. 



II. EFFECTIVE FIELD THEORIES 

In this section, we write the low-energy EFT suited to describe El electromagnetic tran- 
sitions in heavy quarkonia. The EFT is constructed in two steps. The first step consists 
of integrating out modes associated with the scale m. This leads to non-relativist ic QCD 
(NRQCD) coupled to electromagnetism [15|, [l6|. The second step consists of integrating 



out modes associated with the scale mf , which leads to potential non-relativistic QCD (pN- 
RQCp) [l7|, [l8|. The operators of pNRQCD relevant for Ml transitions have been derived 
in [l3|; here we consider the electromagnetic operators responsible for El transitions. 

The EFT matrix elements are counted in powers of f , while the Wilson coefficients of 
NRQCD are series in a^ at the (perturbative) mass scale. Heavy quarkonia may be distin- 
guished in weakly-coupled and strongly-coupled quarkonia. For weakly-coupled quarkonia, 
the binding energy is of the order of or larger than the typical hadronic scale Aqcd- This 
case may be relevant for bottomonium and charmonium ground states. Weakly-coupled 
quarkonia can be treated perturbatively at the typical momentum-transfer scale leading to 
a potential that is Coulombic. The heavy-quark velocity is of order a^ and we may iden- 
tify the momentum-transfer scale with mas and the binding-energy scale with mal. The 
Wilson coefficients of weakly-coupled pNRQCD are series in as-, which can be obtained by 
expanding in as when matching pNRQCD to NRQCD. We will see, however, that, for El 



operators contributing to the transition widths up to relative order f ^, such an expansion is 
not necessary and the obtained results will be valid to all orders in the coupling. In the most 
general weakly-coupled case, the binding energy is not a perturbative scale, ma'^ ~ ^qcd, 
and as at that scale may not be considered an expansion parameter. 

For strongly-coupled quarkonia, Aqcd is larger than the binding energy and possibly of the 
same order as the momentum transfer. This case may be relevant for all higher bottomonium 
and charmonium states. Since El transitions require that at least one involved quarkonium 
state has principal quantum number larger than one, strongly-coupled quarkonia are likely 
to be always involved in such processes. Strongly-coupled quarkonia may in general not be 
treated perturbatively at the typical momentum-transfer scale. This leads, in particular, to 
a potential that is not Coulombic. 

A. NRQCD 

We consider NRQCD coupled to electromagnetism. The Lagrangian has the form 

-^NRQCD = ^2-{ + ^i-f + -Clight • (1) 

The term ^2-f denotes the two-fermion sector of the Lagrangian. If we restrict ourselves to 
operators relevant for El transition widths up to relative order t>^, £2-f reads 

\2 Y\i 
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2m 8m^ 8m^ 

cm sin cm 

2m 8m^ 8m^ 

+ [^ ^ ta'x*, A, ^ -Aj, A^ ^ -A^] , (2) 

where ip is the Pauli spinor field that annihilates a heavy quark of mass m, flavor Q and 
electric charge eeg (e^ = —1/3, Cc = 2/3), and x is the corresponding Pauli spinor that 
creates a heavy antiquark. The gauge fields with superscript "em" are the electromagnetic 
fields, the others are gluon fields, iDo = ido- gT^A'^-eeqA^'^, iB = iV + gT'^A'^ + eegA^"', 
[Dx,E] = DxE-ExD, E^ = F'^, B* = -eijkF^^ /2, W""^ = F^^^"^ and B^^"^ = 
-ei,kF^^^'^/2. 

The coefficients cp, cs, cd, c^, c|™ and c^ are Wilson coefficients of NRQCD. Some 
of them satisfy exact relations dictated by reparameterization (or Poincare) invariance [l9|. 



e.g. 

c7 = 24°^ - 1 , cs = 2cF-l. (3) 

n 

The coefficients are one at leading order, but known at least at one loop [20|]. In particular, 
we have 

4m^l + «:- = i + Cp^ + o{a!), (4) 

where Cp = {N^ - l)/(2iVc) = 4/3 and N^ = 3. The term k"^ is usually identified with 
the anomalous magnetic moment of the heavy quark Q; Kq^ is less than 0.1 for charm and 
bottom. In general, the Wilson coefficients of NRQCD contain also contributions coming 
from virtual photons of energy or momentum of order m. These contributions are suppressed 
by powers of a, the fine structure constant, and will be neglected in the following. 

The terms £4_f and £iight denote respectively the four-fermion sector and the light-field 
sector of the Lagrangian. Light fields include light quarks (assumed to be massless), charm 
quarks in the bottomonium case^, gluons and photons. £4_f and £iight contribute to the 



quarkonium potential and wave functions (see 12| and references therein), but they do not 
provide new couplings of the heavy quarks with the electromagnetic fields relevant for El 
transitions at relative order v"^, which is the accuracy we aim at. 

B. pNRQCD 

In NRQCD, degrees of freedom that scale with the momentum transfer and with the 
binding energy are entangled in physical amplitudes, leading to a non-homogeneous power 
counting. These degrees of freedom are disentangled in pNRQCD, where degrees of freedom 
that scale like mv have been integrated out. The pNRQCD Lagrangian may be decomposed 
in two terms: 

-CpNEQCD = ^0 + 'C7PNRQCD • (5) 

The term Cq denotes the part of the pNRQCD Lagrangian that does not contain heavy- 



^ This is the case when the charm mass is of the order of the momentum transfer. If it is larger, then the 
charm may be integrated out together with the bottom mass, in which case it contributes to the Wilson 
coefficients of NRQCD and does not appear as a light field in the Lagrangian. 



quark couplings to electromagnetisin. In its gauge-invariant form, it reads 

where S = S\c/ y/N], is a quark- ant iquark field that transforms as a singlet under S'f/(3)c 
and ?7(l)em- S is labeled by the quark- ant iquark distance r and depends on the center of 
mass coordinate and time, the derivative V acts on the center of mass coordinate and the 
derivative V^ on the relative distance r. The dots stand for higher-order kinetic energy 
terms. Here and in the following, the trace is meant over color and spin indices. The 
Wilson coefficient Vs-, which is in general a function of r, may be identified with the quark- 
antiquark color-singlet potential. It is organized as an expansion in 1/m, the leading term 
being the static potential, V^ . For weakly-coupled quarkonia, Vs may be calculated in 
perturbation theory, the leading term being the Coulomb potential {Vg = —Cpces/r). For 
strongly-coupled quarkonia, Vs follows from a non-perturbative matching to NRQCD. 

The term A£o describes the propagation of all other low-energy degrees of freedom besides 
the quark-antiquark singlet and their strong interactions. The low-energy degrees of freedom 
depend on the specific quarkonium under scrutiny. For weakly-coupled quarkonia, they are, 
besides the quark-antiquark singlet field, the quark-antiquark field O = v2 0"T", which 
transforms as an octet under SU{3)c and as a singlet under U{l)em, light quarks, low-energy 
gluons and photons; ACq then reads 

ACo = I d'r Tri O^ (iD^ + ^ + Zi: + v<^ O 

J I \ Am m J 

+Va (Qtr ■ (^ES + Str • (7EO) + . . . i + Aight , (7) 

where i^oO = idoO - ^[T'^Ag,0] and iDO = iVO + ^[T-^A^O]. Like S, the fields C^ 
are labeled by the quark-antiquark distance r and depend on the center of mass coordinate 
and time. To ensure that gluons and photons are of low-energy (i.e. carry energy and 
momentum lower than the typical momentum transfer in the quark-antiquark system) gluon 
and photon fields are multipole expanded in the relative distance r and depend only on the 
center of mass coordinate and time. Hence the Lagrangian is organized as an expansion 
in r and l/m (inherited from NRQCD). The dots in (JTj) stand for terms that contribute 
to El transitions beyond our accuracy. The Wilson coefficient Vq niay be identified with 
the quark-antiquark color-octet potential. It is organized as an expansion in l/m, the 



leading term being the static potential, Vq . At leading order in perturbation theory, 



Vq = as/{2Ncr). The Wilson coefficient Va = 1 + 0{a'^) [2l| provides the strength of 
the chromoelectric dipole interaction. £iight denotes the light-field sector of the Lagrangian; 
light fields include light quarks (assumed to be massless), gluons and photons. For strongly- 
coupled quarkonia, after having integrated out Aqcd, only degrees of freedom that are color 



singlet are possible |22|, |23|. These are, besides the quark-antiquark color-singlet field and 
photons, the Goldstone bosons associated to the spontaneous breaking of chiral symmetry. 
The effects of Goldstone bosons on El transitions go beyond our accuracy and will be 
neglected. Hence, for strongly-coupled quarkonia we set A£o = —F'^F^'^'^^ /A, whereas all 
the complication of the non-perturbative treatment goes in the determination of V5. 

The term /Z-ypNRQCD describes the coupling of heavy quark-antiquark pairs with low- 
energy photons, like those responsible for electromagnetic transitions. The power counting 
goes as follows 

r ~ l/mv, Vr ^mv, V ~ mv'^ , k^, W^, W^ ~ A;^ ; (8) 

in the case of weakly-coupled quarkonia, one has to consider also low-energy gluons that 
scale with mv'^ or k^. The leading operator responsible for El transitions is the electric 
dipole operator S^r ■ eegE^^S, while operators relevant for El transitions at relative order 
f^ are those suppressed by t>^ with respect to it. The part of £^pNR Qcn relevant for El 



transitions is (operators relevant for Ml transitions have been listed in 13l |) 



-C^pNRQCD = j dh TrJ r^-^S^r ■ eeQE^"^S 



+ ^\/("^)''-^SV ■ [(r • V)2eeQE^'"]S 

+ ^V^v..(rx{rV)B)gt |v^.^ r X [(r ■ V)eeQB'^^]}S 
12m 

Am 

_yHE^v,.)p^ 0-] ■ (eegE-^ X V,)S \ + A£^pnrqcd • (9) 



Am? 

Note that the condition k^ <^ mv guarantees that we can multipole expand the electro- 
magnetic fields regardless of the weakly- or strongly-coupled nature of the quarkonia. On 



symmetry grounds, more terms than those hsted in ([9j) are possible. However, as we will 
argue in the next section, these are the only ones that get contributions from matching 
with NRQCD. The first line contains the leading electric dipole operator, all other operators 
are suppressed by w^. The coefficients V''^ , y(rVfr-E ^ yv-(rxB)^ •^v,{rx{rV)B)^ y{rV)a-B ^^^ 
Ya-{ExVr) QJ.Q Wilson coefficients that will be computed in the next section. 

The term A£^pnrqcd contains the electromagnetic couplings with the other low-energy 
degrees of freedom besides the quark- antiquark singlet. For weakly-coupled quarkonia the 
only relevant term for El transitions at relative order f ^ is 

A£^pNRQCD = j d\ Tr [VS^O^v ■ eeQE^^^O} , (10) 

which is the electric dipole operator for quark-antiquark states in a color-octet configuration. 
The corresponding Wilson coefficient is Vq^ . For strongly-coupled quarkonia, we may set 

A£^pNRQCD = 0. 

III. MATCHING 

We calculate here the Wilson coefficients of -C^pnrqcd, while the Wilson coefficients of Cq 
can be found in [1^ and references therein. Calculating the Wilson coefficients of -C^pnrqcd 
requires to match NRQCD quark-antiquark Green's functions coupled to one external elec- 
tromagnetic field with pNRQCD ones. The electromagnetic field can be multipole expanded. 



The matching can be done order by order in 1/m |12l . |20|. We will perform the matching 
at leading order in the electromagnetic coupling and at all orders in the strong coupling. 
Suitable Green's functions are static Wilson loops with electromagnetic and gluon field in- 
sertions. 

A. Matching at 0{al) 

Before calculating the Wilson coefficients of -C^pnrqcd to all orders, we consider, in the 
case of weakly-coupled quarkonium, the matching from NRQCD to pNRQCD at 0{a^). 
The calculation can be performed by expanding and redefining the fields in the NRQCD 
Lagrangian. This goes in two steps. 



(i) First, one projects NRQCD on the quark-antiquark Fock space spanned by 

/ d^xi(fx2ij\xi,t)(p{xi,X2,t)x{x2,t)\\ight) , (11) 

where f{xi, X2, t) is a 3 ® 3 tensor in color space and a 2 ® 2 tensor in spin space, and 
|hght) is a state that contains an arbitrary number of low-energy gluons, photons and 
light quarks, but no heavy quarks. 

(ii) Second, one decomposes 

(p{xi,X2,t) = Fexpiig / A ■ rfx J S'(r,R, t) 

+Pexplig A-dx] 0'{r,R,t) explig A-rfxJ, (12) 

S'(r, R, t) = exp (ieeg f A""' ■ rfx j S(r, R, t) , (13) 

0'(r, R, t) = exp (ieeq f A^'" ■ rfx j 0(r, R, t) , (14) 

where P stands for path ordering, R = (xi + X2)/2 is the center of mass coordinate, 
and r = Xi — X2. The decomposition ensures the gauge invariance of the pNRQCD 
operators. Finally, all electromagnetic fields are multipole expanded in r. 

The resulting 0{a^) expression of the Wilson coefficients is^ 

yr-E _ yr-E _ y{r\')'^r-E _ y\'-{rxB) _ T/Vr-(rx(rV)B) _ i /i r\ 

y{r^)a-B ^ ^cm^ ^^g^ 

Surprisingly, these relations will turn out to be vahd to all orders. 

B. Matching photons coupled to light quarks 

Photons may couple to heavy quarks or to light quarks. If we treat the u, d and s quarks 
as massless, then the QCD Lagrangian is S'[/(3)-flavor symmetric and the three quarks only 



^ 0{a^) refers to the matching between NRQCD and pNRQCD: the Wilson coefficients of NRQCD are 
kept unexpanded. 




FIG. 1: Leading-order diagram with a c-quark loop coupled to an external photon in NRQCD. 
Dashed lines stand for longitudinal gluons. 

differ in the electric charges. Since the sum of the three charges vanishes, so does the sum 
of all diagrams where the photon couples to the three massless quarks. 

In the bottomonium case, if the charm-quark mass is of the same order as the typical 
momentum transfer in the system then it should be integrated out with that scale. Indeed 
it contributes to the potential |2J], although in the spectrum it appears to decouple [25|. 
The leading-order diagram with a c-quark loop coupled to an external photon in NRQCD 
is shown in Fig. [TJ It is of order a^{mc) ~ v^, thus beyond our accuracy. 

C. Matching of reducible diagrams 

It is useful to observe that certain classes of diagrams in NRQCD contribute just to 
reducible diagrams in pNRQCD, i.e. tree-level diagrams made of potentials and electro- 
magnetic operators of lower order. Hence, they do not contribute to the matching of new 
operators. This happens whenever a gluonic or electromagnetic contribution can be factori- 
zed out. In the following, we will illustrate some of these cases. 

(i) Consider static NRQCD. It describes the propagation of a static quark located in xi 
and a static antiquark located in X2. We call \n), states that contain an arbitrary 
number of low-energy gluons and light quarks, but no heavy quarks. Furthermore we 
assume that '?/'^(xi)x(x2)|n) is an eigenstate of the static NRQCD Hamiltonian with 
eigenvalue En- The eigenstates are normalized such that {n\n) = 1. In general, both 
\n) and En depend on r. 
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(a) The matching between NRQCD and pNRQCD follows from equating in the large time 
limit the time-ordered NRQCD amplitude 

(0|^(xi)x^(x2) /"rfSi {^\z^)OM)^{z,) + c.c.) ■ ■ • 

X jd'z^ {^\z,,)On{ZnMZn) + CC.) ^^ (x;);^(x^) |0) 
= (0| /"rfti(Oi(ti)+C.C.)--- fdtn{On{tn)+ c.c.) 

xOih - ts) ■ ■ ■ 0{tn-i - t„)|0) <53(xi - x;) 6%^2 - x'2), 
with the corresponding pNRQCD one. The NRQCD amplitude may be identified in the 



large time limit with a static Wilson loop with Oi, ..., 0„ and c.c. field insertions 2^, 



231]. The operators On{t) are gluonic and electromagnetic operators in the Heisenberg 
representation, c.c. stands for charge conjugation (after the equality, it stands for 
charge conjugation and Xi f-)- X2 exchange). The corresponding pNRQCD amplitude 
follows from '?/'''"(xi)x(x2)|0) — )■ 5'"''(r, R)|vac) and Eq = Va , where |vac) is the vacuum 



22 



of pNRQCD normalized such that (vac | vac) = 1, see 

(Hi) In particular, suppose that the non-electromagnetic operator tp'^Oip + c.c. of NRQCD 
matches the operator / d'^r" S'^5VS of pNRQCD, then the matching condition reads 

(0| f dt (0(X1, t) + 0^(X2, t)) |0) 53(X1 - X'l) 5=^(X2 - x'2) 

= (vac|5(r,R) j dt f d^R" f d^r" S^5VS ^^(r', R')|vac), 
where O^ is the charge conjugated of O. This implies 

fdte'^"' (O|0(xi) + O^(x2)|0)e-^^o* = f dte'^s°^' 6V e-'^s°^\ 

and eventually 

(0|O(xi) + O^(x2)|0) = 5l^. (18) 

Another example is the case of a NRQCD amplitude with insertions of two non- 
electromagnetic operators ip'^Oiip and ip'^02ip with their charge-conjugated partners. 
Such an amplitude matches the pNRQCD reducible amplitude made of two insertions 
of 51^1 = (0|Oi+O^|0) and 51^2 = (0|O2+Oi|0), and the pNRQCD amplitude associated 
with a new potential 5Vi2, 

/■CXI 

5Vi2 = -M dt ^(OlOi + 0{\n)e-'^-\n\02 + O^|0) 6*^°* + (1^2); (19) 



n^O 
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n 7^ intermediate state contributions that are exponentially suppressed at large times 
have been neglected. 

(iv) Consider now an electromagnetic operator of NRQCD, ip'^O'^'^ilj + c.c, that matches 
the operator / d^r" S^V'^S of pNRQCD. If O"^ commutes with the gluon fields (at 
least to the order in the power counting we are interested in), then O*'™ commutes 
with the static NRQCD Hamiltonian and its eigenstates (e.g. [O^'^, |0)] = 0). The 
argument of paragraph (Hi) then implies 

Qcm ^ Qcmc ^ jycm _ (20) 

(v) A simple extension of the previous case is the one of a NRQCD amplitude with Oi, 
..., On, O^™ field insertions under the condition that O*'™ commutes with all gluonic 
operators Oi, ..., On (at the order in the power counting we are interested in). The 
amplitude is then proportional to 



Ol(ti)02(t2) . . . 0„(t„)0""(tem) [0{tcrn " tl)^(ti - ta) • • • ^(4-1 - ^n) 

+6[ti — tcm)6[tcm — 12) ■ ■ ■ 0{tn-l — t„) + . . . 
+9{ti —t2)... 6{tn-l — tn)0{tn — tcm) j 

= Oi(ti)02(t2) . . . On{tn)e{h - ^2) . . . ^(t„_i - t„) 0^'"(te„,), (21) 

where the equality follows from the recursive use of 9 {tcm— 1')0 (t' —t") +9 {f —tcni)0 {tem — 
t") = 9{t' — t")9{tcni — t"). The equality implies that the amplitude reduces to the 
product of a pure gluonic amplitude and the electromagnetic vertex O'^™. Therefore 
it matches the reducible diagram of pNRQCD made of some potentials (those that 
match the non-electromagnetic NRQCD amplitude with Oi, ..., 0„ field insertions) 
and the electromagnetic vertex p^™. 

(vi) We match now the following amplitude of NRQCD: 



(0| / dti / dt2 (0(xi,ti)+C.c)(0'="^(xi,t2)+C.c)^(tl-t2) 

+ (0^-(xi, h) + c.c) (0(xi, t2) + c.c) 9{h - t2) |0) (5-^(xi - x'^) 6%^2 - x^), (22) 

under the general assumption that O'^™ does not commute with O. By inserting a 
complete set of eigenstates and making use of the fact that 0°™ commutes with the 

12 
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FIG. 2: Matching of the amplitude (vi): before the arrow are the NRQCD diagrams, after the arrow 
the pNRQCD ones. The continuous hne in the pNRQCD diagrams stands for the quark-antiquark 
singlet propagator. 



gluon fields, see paragraph (iv), the amplitude becomes 
fdh /"dt2^(ti -t2)e'^«*^| (O|0(xi) + O^(x2)|0)e-'^«(*^-*^) (0^'"(xi) + ©^'"^(xs 

+ (0^'"(xi) + 0^"^^(x2)) e-'^«(*^-*^)(0|O(xi) + O^(x2)|0) 

xe-^^«*^53(xi-x;)53(x2-x'2), 

which, according to (fT8|) and (120|) . matches the reducible pNRQCD amplitude 



)) 



f dh f dt2 9{h - ts) e'^i°'*i I 6V e-^^i°^(*i-*2)p«'- + p-- e" 



■iV^°\h-t2)^y 



xe-*^s '*2 53(xi - x'J 6%^2 - x'2) . 

Graphically this is shown in Fig. [2J Therefore an amplitude like Eq. (122 p cancels in 
the matching against a reducible amplitude of pNRQCD. 

(vii) Finally, let us consider the case of a NRQCD amplitude with two gluonic operators 
Oi and O2 and corresponding charge-conjugated operators, and an electromagnetic 
operator 0°™ that commutes with one of the gluonic operators, say O2. The amplitude 
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reads 



f dh I dt2 /"rft3e(ti-t2)e(ti-t3)e*^o*i 

xe-'^o*2<53(xi-x;)53(x2-x'2) 
x$^ J (O^-(xi) + 0'='"^(x2)) (O|0i(xi) + OJ(x2)|n)e-^^"(*^-*^)(ri|02(xi) + O2^(x2)|0)! 



xe-'^"*2 53(xi-x;)53(x2-x'2) 



+ (1^2). (23) 

The n = component of the amphtude matches the reducible pNRQCD diagrams 
made of insertions of the potentials 5Vi, 6V2 and of the electromagnetic vertex 'P*'™. 
The sum of the n 7^ components matches 

f dti f dt2 e{ti - t2) e^^s"'*! I iSVu e-'^S°\il-i2)'pem ^ pem ^-ivf it,-t,)-^y^^ I 



_,-T/(0). 



Therefore, also this kind of amplitude does not induce new operators in pNRQCD. 

D. Matching at 0(1/^°) 

Amplitudes that contribute to the 0{l/m^) matching of £7pnrqcd may contain the 
NRQCD operators (see Eq. ^) 

-^^eegAl^^, -ij^gAo^, (24) 

and the corresponding c.c. ones. Although they just contain one electromagnetic operator, 
they may contain an arbitrary number of longitudinal gluons for they are not suppressed by 
any power of 1/m. The electromagnetic field Aq™ commutes with the gluon fields. Therefore 
it satisfies the condition of Sec. IIIICl paragraph (iv), and the matching condition is given by 
Eq. ( |20l) . It tells that electromagnetic operators of order l/m° do not get QCD corrections. 
Hence all C(l/m°) Wilson coefficients of -C^pnrqcd are fixed at their C(as) value: 

yr-E ^ y{rV)^vE ^ ^ _ ^25) 

14 



For weakly-coupled quarkonia, also the quark-antiquark color-octet sector is relevant. In 
this case, the matching is performed by considering NRQCD amplitudes between initial and 



final states that match the color octet state in pNRQCD 18|. Since the electromagnetic 



field commutes also with these states, the result of the matching to all orders is 

V^-^ = 1 . (26) 

Finally, we observe that, while our argument fixes V^'^ to all orders, the same argument 
does not apply to the Wilson coefficient Va of the chromoelectric dipole operator in Eq. ([7j). 
The reason is that the field Aq does not commute with the gluon fields. 

E. Matching at 0{l/m) 

Amplitudes that contribute to the 0{l/m) matching of /^-ypNRQCD contain, besides an 
arbitrary number of operators of order l/m^, one of the NRQCD operators of order 1/m 
(see Eq. (ED): 

■rj2 em 

^t ^, -L^V.^B^, -^^V . eegB^-^, (27) 

2m 2m 2m 

and the corresponding c.c. ones. We call the first operator the kinetic energy operator, 
the second one the chromomagnetic dipole operator and the third one the magnetic dipole 
operator. The 0{l/m) amplitudes fall in one of the following categories. 

1. The photon is coupled to the magnetic dipole operator. This kind of diagrams matches 
spin-dependent operators and may contribute to \/(^''^)<^-S_ Using the same argument as for 
the 1/m^ matching, the amplitude factorizes. Therefore, it holds to all orders that 

y{rV)a-B ^ ^em ^ ^28) 

and no other operator is generated. 

2. The photon is coupled to the kinetic energy operator (see e.g. Fig. [3]). This kind 
of diagrams matches only operators with magnetic fields, hence they may contribute to 
yV(rxB) ^^^ ^v,.(rx(rV)B)^ Siucc the electromagnetic coupling is embedded in a covariant 
derivative operator, 

V^f(V - ieeQA^'^)V^ + c.c. , (29) 
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FIG. 3: Examples of diagrams at order 1/m, where the electromagnetic coupling is contained in a 
covariant derivative. 

this implies that diagrams involving one kinetic energy operator match pNRQCD operators 
of the form^ 

Jd'r S'^ ('|^,(V.^ -^eeQA^-(xi,t))4 +C.C.') S' 

+ fd'rS'^ (!.^-^r■,l{V,,-^eeQA'^{x^,t))\+c.c^S'. (30) 



After the field redefinition flT3|) and having multipole expanded the electromagnetic fields 
up to include order r^ contributions, the pNRQCD operators may be rewritten as 

+ L-VS* |^r-,2'iV,|s. (31) 

Switching off the electromagnetic interaction, this expression should match Eq. (J6]), which 
fixes 6Va = 1/4 and 6Vb = to all orders in perturbation theory. The fact that the kinetic 
energy in pNRQCD is protected against quantum corrections is a direct consequence of 
Poincare invariance |26|]. We conclude, therefore, that 

yV-irxB) ^ yVr-(rx(rV)B) ^ ^ /gg-) 

holds to all orders in perturbation theory. 

3. The photon is longitudinal and coupled to the heavy quark/antiquark lines. Such 
diagrams match only operators with electric fields. If we consider diagrams with one insertion 



^ For the purpose of the present discussion, we do not consider operators that do not depend neither on 
the momentum nor on the electromagnetic field. They contribute to the 1/m potential and have been 
analyzed in ;22J . 
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FIG. 4: Examples of diagrams at order Xjm with an external longitudinal photon. 

of a chromomagnetic dipole operator or a kinetic energy operator (for the latter see Fig. H]), 

then such diagrams may possibly contribute at order \jm to operators that vanish at (9(ag), 

for instance, 

[ d^r —S^r-eeoB'^'^S. (33) 

J mr 

In both cases, the argument developed in Sec. IIIICl paragraph (vi), applies, for Aq™ com- 
mutes with gluons. Hence, all such diagrams cancel in the matching against reducible 
pNRQCD diagrams, and an operator like the one written above does not show up even at 
higher orders in the coupling constant. 

F. Matching at 0{l/m^) 

Diagrams that contribute to the matching at order l/m? contain, besides an arbitrary 
number of operators of order l/m°, either two operators of order 1/m or one of the following 
operators of order 1/m?: 

cm cm 

which we call the (chromo)electric spin-orbit operator and the (chromo)electric Darwin 
operator respectively. Since many different diagrams can contribute, it is convenient to 
distinguish between spin-dependent and spin-independent amplitudes. It is always implicitly 
assumed that diagrams may contain an arbitrary number of l/m° gluonic vertices. 

1. Spin- dependent diagrams 

First, we consider the matching of spin- dependent operators. 
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FIG. 5: Some diagrams that contribute to spin-dependent operators of order Xjrr? containing an 
electric field. 

1. We consider diagrams made of two chromomagnetic dipole operators and a longitudinal 
photon. Since the photon commutes with the chromomagnetic dipole operators, the result 
of Sec. unci paragraph (v)^ applies: these diagrams do not contribute to the matching of 
new operators. 

2. Diagrams that contain one magnetic dipole operator and a chromomagnetic dipole 
one are of the type discussed in Sec. IIIICl paragraph (vi). Hence they do not contribute to 
the matching of new operators. Moreover, in [13| , it has been pointed out that this kind of 
amplitudes, proportional to the expectation value of a chromomagnetic field, vanishes for 
parity. 

3. The situation is similar for diagrams that contain one magnetic dipole operator and 
a kinetic energy one. Since the magnetic dipole operator commutes with the gluons, we are 
in the situation of Sec. IIIICl paragraph (vi)^ and this type of diagrams does not contribute 
to the matching of new operators. 

4. Also diagrams that contain one chromoelectric spin-orbit operator and one longitu- 
dinal photon (see e.g. diagram 1 in Fig. [5]) fall under the situation discussed in Sec. IlllCt 
paragraph (vi)^ and do not contribute to the matching of new operators. 

5. We consider now diagrams made with one chromomagnetic dipole operator, one kinetic 
energy operator and a longitudinal photon (see e.g. diagrams 2 and 3 in Fig. |5]). Since the 
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longitudinal photon commutes with the chromomagnetic dipole operator, the argument of 
Sec. IlllCt paragraph (vii), applies and therefore such diagrams do not contribute to the 
matching of new operators in pNRQCD. 

6. Diagrams that contain one electric spin-orbit operator (see e.g. diagrams 4 and 5 in 
Fig. E]) contribute to operators that are at least v"^ suppressed with respect to the leading 
electric dipole operator of pNRQCD. These operators depend on spin and on an electric 
field. For time inversion invariance, they must also depend on the quarkonium momentum. 
Hence, only the derivative part of the covariant derivative of the electric spin-orbit operator 

pern 

contributes and only if it does not act on the gluon fields. Since —^cr ■ [— 'iVx,eeQE^™l 
effectively acts as an operator that commutes with the gluons, we are in the situation 
discussed in Sec. IIIICl paragraph (iv), that led to the matching condition ( l20i) . In our case 
and at leading order in the multipole expansion, the matching condition becomes 

which is exact. Operators that come from higher-orders in the multipole expansion do not 
need to be considered here because they are beyond our accuracy (although the matching 
fixes also the Wilson coefficients of those operators to all orders). Operators that are possible 
for symmetry arguments alone, like for instance 

1 



^2^2 



[S\cTmrxtVry,r-eeQE'^}S, 



cannot be generated at any order in the strong-coupling constant and, therefore, may be set 
to zero in pNRQCD. Finally, we observe that the observation made in 22| that (0|D^.J0) = 
Vxi and (0|Da;2|0) = V^jj would lead to the same result. 

7. Diagrams containing either one chromomagnetic operator and a kinetic energy op- 
erator coupled to an external photon or one chromoelectric spin-orbit operator with the 
electromagnetic field encoded in the covariant derivative have been calculated to all orders 



in 131] with an argument similar to the one used in Sec. IIIIEt paragraph 2. They contribute 



to the Wilson coefficient of one single operator. 



^ 1 {St,^.[rx(rxeeQBn]}S, (35) 



which, however, is relevant only for Ml transitions and not for El ones, since it does not 
change the parity of the quarkonium state. 
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FIG. 6: Some diagrams contributing to the matching of spin-independent operators of order \jrv? . 
2. Spin-independent diagrams 

Now we consider the matching of spin-independent operators. 

1. Diagrams with one insertion of an electric Darwin operator do not contribute beyond 
0{a^) to the matching, for the electric Darwin operator commutes with gluons and the 
conclusion of Sec. IIIIC[ paragraph (iv), applies. These diagrams match into the electric 
Darwin operator of pNRQCD. Such an operator has not been displayed in Eq. (|9]), because 
it does not contribute to El transitions. 

2. Diagrams that contain one chromoelectric Darwin operator and a longitudinal photon 
are of the type discussed in Sec. IIIICl paragraph (v). The electromagnetic interaction 
factorizes and the contribution cancels in the matching. 

3. We consider diagrams containing two kinetic energy operators. The (transverse) 
electromagnetic field is embedded in one of the covariant derivatives. These diagrams match 
spin-independent operators of pNRQCD with one magnetic field. Because of the pNRQCD 
symmetry under charge conjugation and Xi ^ X2 exchange, the allowed spin-independent 
operators with one magnetic field must contain an odd number of r or V^. Moreover, 
because of parity invariance, at least one center of mass derivative has to be present too. 
However, such operators are too much suppressed to be relevant for El transitions at relative 
order f^, an example being the operator 

1 



m?r 



SHzV^rxeeQB^'^lS 



which is of relative order v^ with respect to the leading electric dipole operator. 

4. Finally, we consider diagrams containing two kinetic energy operators and an external 
longitudinal photon, see Fig. El We observe that A^™(xi) + ^^^^(xa) = r • VAg"" + 0{r'^). 
Hence, these diagrams contribute to operators that are at least l/rn? x {mvY ~ v"^ suppressed 
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with respect to the leading electric dipole operator of pNRQCD; operators generated by 
higher orders in the multipole expansion are beyond our accuracy. Let us consider now the 
kinetic energy operators. If one of the derivatives acts on the longitudinal photon then the 
diagrams contribute to operators that are at least v^ suppressed with respect to the leading 
electric dipole operator of pNRQCD. Operators of relative order v^ are beyond our accuracy. 
If one of the derivatives acts on the factor r in r ■ VAq™ then the other derivatives do not. 
As a consequence, the electromagnetic operator r ■ VAq™ commutes at relative order v"^ 
with at least one of the kinetic energy operators and we are in the situation described in 
Sec. IIIIC[ paragraph (vii). We conclude that in pNRQCD these diagrams cancel against 
iterations of lower-order potentials and electromagnetic vertices and do not contribute to 
any new operator. In particular, possible operators allowed by the symmetries of pNRQCD, 
like, for instance, 

cannot be generated at any order in the strong-coupling constant. 

G. Concluding remarks 

In this section, we have matched non-perturbatively all operators of pNRQCD relevant to 
describe El transitions at relative order v"^. It turns out that the situation for El transitions 
is different from the one for Ml transitions first discussed in |lc 

For El transitions no 0{l/im?) operator is relevant at relative order f ^, as it is shown in 
Eq. (I9]). In particular, we can neglect the matching of NRQCD diagrams with three kinetic 
energy operator insertions or with one insertion of the operator D^/(8?7i^). 

For Ml transitions between strongly-coupled states, corrections of relative order v"^ require 
the non-perturbative matching of 1/iv? operators. At this order one can perform an exact 
matching for all but one relevant operator (see appendix). This is the operator 

1 



rn?r'^ 



V^ {S^o-■eeQB^'"}S, (36) 



whose Wilson coefficient, V^ , can possibly get QCD corrections at the momentum-transfer 
scale. These are encoded in the expectation value of some suitable Wilson loop, whose 
explicit expression is at present unknown. 
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IV. El TRANSITIONS 

After having derived the relevant pNRQCD Lagrangian, we will now proceed in the 
calculation of the electric dipole transition rates, starting from the non-relativistic limit. This 



will be used to fix the notation and discuss the wave functions. We will follow closely |13|. 



A. Quarkonium states 

We consider radiative transitions, if — )■ if'7, between a quarkonium H and a quarkonium 
H' . A quarkonium state |-ff(P, A)) is an eigenstate of the pNRQCD Hamiltonian with the 
quantum numbers of a quarkonium H with polarization A. We normalize it in the non- 
relativistic way, i.e. 

{H{V\\')\H{V,\)) = (27r)3<53(P-P')(5,v. (37) 

The leading-order quarkonium state is defined as 

|i7(P,A))(°) = I d?R /"rfVe'P-^Tr|</)gJ^)(r)S^(r,R)|vac)} , (38) 



which is also an eigenstate of the total spin, the orbital angular momentum and the center 
of mass momentum P of the quark- ant iquark pair. The wave function (|y^\)^ (r) is an eigen- 
function of the leading pNRQCD singlet Hamiltonian /i^ , i.e. a solution of the Schrodinger 



equation 

,(0)^(0) _ / _ V^ ^ ,^(0) ^ ^(0) _ p(0)^(0) 

m 
The eigenvalue E)j is the leading-order binding energy. 



''™<!., ^[~ + vr] <U = ^«'<°!., . (39) 



The wave functions have the following angular structures for L = 13| 







3^^no(r), (40) 



nkwi^) = J^ Rnoir) o- ■ e„35,(A) , (41) 
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for L = 1 



.(0) .r 




niPi(A)VV Yg^ 

3 

n3Pi(A)iM - YTevr 



i?„i(r)e„ip,(A)-f, (42) 



^°\/r) = ^/-i?„i(r)a-r, (43) 

'^"^ W = \/T^i?ni(r)o--(rxe„3p,(A)), (44) 



^it.(A)(r) = V 1^ i?ni(r) cr^h^^{X)r^ , (45) 



and for L = 2 y, y 



,(0) . . / 15 



n^^.(A)W = ^Y^Rn2{r)rX.^^{X)r\ (46) 



^0) 



1 



n3D,(A)(r) = V Te^^'^^^'') [3(e„3,,,(A) ■ r)(^ ■ f ) - ^ ■ e„3^,(A)] , (47) 




'n^DiWy-J — V g^ -'^"^V ;^ "-nazja 



^"^ ^r) = \li-Rn2ir) rX^^Sm X ay , (48) 



<^%.w(r) = ^^ i?„2(r) FFif:?^3(A)^^ . (49) 

The vectors e„25+ii_,(A) denote orthonormal polarization vectors of the quarkonium state. 
The tensors h^^2s+iL (A) and -ff^i-s+i^ (A) are completely symmetric, traceless (the tensor 
H^^^ has vanishing partial traces, i.e. H^2s+il = 0) and normalized as 

h^.^L,ims.^LA^) = ^AA' = //:^:.,^(A)iC^,,,^(A') . (50) 



Whereas the quarkonium state is normalized in the non-relativistic fashion f l37j) . the one 
photon state, |7(k, cr)), is normalized in the usual Lorentz-invariant way 

{^{k,a)\^ik',a')) = 2k{2nf5%k-k')5^^,. (51) 

This implies that external electric or magnetic fields project on a one photon state as 

(7(k,(T)|E'='"(R)|vac) = -ike*{a)e-'^-^, (52) 

(7(k,(T)|B^'°(R)|vac) = -ikxe*{a)e-'^'^, (53) 

where e(o") is the photon polarization vector (the dependence of e(cr) on k is understood). 
The photon transversality requires e(a) ■ k = 0. 
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Ph = {Mh, 0) X. P^. = U^ + M-fj„ -k) 



FIG. 7: Kinematics of the transition H ^ H' 'j in the center of mass frame. 
Finally, the quarkoniuin and photon polarizations satisfy the relations 

E<-+^A.(^)<-+^L,(^) = ^'^ (54) 

A 

J2 K:s,.:^^{\)h%,.^^{\) = ^i5V + 5V) - ^S^^'^ , (55) 

A 

+ permutations of the indices i, j, k and l,m,n] , (56) 
J2 e'(ff)e^'*(ff) = S'^ - k*k^' . (57) 

B. Transition amplitudes and rates 

In the rest frame of the initial quarkonium, see Fig. [71 the transition amplitude from a 
quarkonium H with polarization A to a quarkonium H' with momentum P' and polarization 
A', and a photon with energy 

*' - W - 4m^ - ^''« - "'"^ ^ Kll) ■ <'*' 

and polarization a is given by 

^ [H{0, A) ^ //'(-k, A')7(k, (t)] {2nf5%P' + k) = 

-{H'{P',X')-f{k,a)\jd^RC,p^nQCB\HiO,X)). (59) 
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The corresponding transition width reads 

<PP' d^k 1 



-"- - 1 ^^wr^Tk p')'***^" - ' - ^') w. 5:m«.-. 



1 / k 

'- ' "'^" ^ ■' (60) 



^) Hdhk fdnik)5ik-k^)^y2\AH^H' 

IhJ Jo J ^Af-r^ 



AA'o- 

where the initial state is averaged over the polarizations, whose number is A^a, and Au^H'-i 
is a short-hand notation for A [H{Q^ A) — ?■ i7'(— k, A')7(k, a)]. 

C. The non-relativistic limit 

The leading operator responsible for El transitions in pNRQCD is the electric dipole 
operator 

Cei = I dh Tr {SV ■ eegE^'^S} . (61) 

From Eqs. fl38l) . fl52|) . fl59l) and the equal-time commutation rules for singlet fields (spin 
indices are written down explicitly) 

4.(r, R, t), 5,z(r', R', t)] = (5,z5,fc(53(r - r')(5^(R - R'), (62) 

it follows that the leading amplitudes for the El transitions n^Pj — )■ n''^S'i7 and n^Pi — >■ 
n'^So'j read 

^SVo^.'35.. = ^/3(nl -> n'O) e*{a) • e„,35,(A0 , (63) 



■^%,^n'^s,, = -^h{nl ^ n'O) e*{a) ■ (e„3p,(A) x e„,35,(A')) , (64) 



,(0) _ ieeqk 

^Sp.^n'35,, = ^^3(nl ^ n'O) e:,,3,,(A')/.:^3p,(A)6*^(a) 



(65) 

'^r/p.^n-so. = ^^3(nl -> n'O) 6*(a) • e„.p,(A) , (66) 

where we have defined 

POO 

lM{nL^n'L')= / rfrr^ /?„a'(r) i?„L(r) , (67) 

r°° d 

I'j^inL^n'L') = / drr''Rn>L'ir)—RnLir). (68) 
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The second definition is for furtlier use. We also give tlie amplitudes for n^Di — )■ n'^Pj'j: 
•^nk^n'sPo, = '-^^^ ^3(r.2 ^ n'l) e*(a) • e„3^,(A) , (69) 

^Sk-n'3p,, = ^^3(^2 -> n'l) e*{a) ■ (e„3^,(A) x e„,3p,(A')) , (70) 

^X^n'^P.. = ^^3(^2 ^ n'l) K,nMn^pS>^')e*K^) ■ (71) 

From (1601) and the relations ( I54ll571) . it follows that the non-relativistic decay rates are 



r(o) _ p(o) A _ ^7 ^ /72^ 

p(0) _ p{0) A _ ^7 ^ /yoN 

with 

r™-5a4*?|/3(nl^n'0)P~^. (74) 

El transition rates are of order k^/{m^v'^), which means that they happen 1/f^ more fre- 
quently than allowed Ml transitions at the same photon energy. The terms proportional to 
k^/M^spj ~ k^/Mnip^ ~ mv"^ /{2m) are suppressed by v"^ and can therefore be neglected at 
leading order. They contribute and have to be accounted for at relative order f ^, which we 
will compute in the next section. We note that at leading order the transition rate is inde- 
pendent of the initial total angular momentum J . The reverse transitions depend instead 
on the final total angular momentum: 

.(0) _2J + 1 (0) / k^ 



^n^s.^w'^Pj,- 3 r„„, (^1 ^^^^J > (75) 

These results agree, neglecting k^/m corrections, with the general non-relativistic for- 



mula |28| 



2 



T%^,^^^^as^,^,^^^ = \ae%k'^[h{nL ^ n'L')f {2J' + l)ms.^mL') | ^, ^ "[ ^ (77) 



'J' 



where the last term in the brackets denotes a Wigner 6-j coefficient. 

V. RELATIVISTIC CORRECTIONS TO El TRANSITIONS 

One of the main advantages of using pNRQCD is that it accounts for the corrections to 
the decay amplitudes in a systematic fashion. We will split the analysis into a part that 
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accounts for the electromagnetic interaction terms in the pNRQCD Lagrangian suppressed 
by (9(f^) with respect to the leading electric dipole operator fl6T]) . and a part that accounts 
for f ^ corrections to the quarkonium state. We will concentrate on the decay rr'Pj — )■ n'^Si 7 
to compare with the result in [id]. The extension to other processes like n^Pi — t- r2'^S'o7, 
rr'Si — !■ n'^Pj'j and u^Sq — )■ n'^Pi'y is straightforward. Final results for all these radiative 
transitions will be given in Sec. I VII For a non- vanishing leading-order transition amplitude 

•^n2s+iLj(A)^n'2s+iL'j,(A')7' ^^ define 

A. Corrections induced by El operators of relative order v'^ 

At subleading order in the decay rate, i.e. at order k^/m?, all the interaction terms 
displayed in the Lagrangian ([9]) beyond the leading electric dipole operator contribute. 



(1) The correction induced by the operator S'^r ■ [(r ■ V)^eeQE''™]S'/24 is 

Id) _ fc^ h{nl ^ n-0) 

•^n3p,^n'35i7 120/3(^1 -^ n'O) • ^ ' 

(2) The correction induced by the operator iS\^-,v x eeqG'^^} S / {Am) is 

(3) The operator iS'''"{ Vr-, r x [(r • V)eeQB''™]}5'/(12m) corrects the leading-order ampli- 
tude by an amount 

(4) The spin-dependent operator c^[S'''', cr] ■ [(r ■ V)eeQB°™]S'/(4m) contributes with 



„em 



- , K'r' — 

A3Po^n'35i7 = -1^ ' (82) 

^(4) ^ fccf° (k ■ e„,35je* ■ (k x e„3pj 

'^n3P,^n'35,7 2m 6* ■ (e„,3^, X e„3pj ' ^^^^ 



A. 



(4) _ '^'-p 



[e*-e„,35jk*-(k-e„,35je* 



/i%pk^' 

n^P2 



„3p,^„,35,^ 2m e*'h%^e\,^ 
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(84) 



where the quarkonium and photon polarization numbers are not exphcitly shown. 
Using Eqs. f l54|) -f l57|) the corresponding decay rates yield 

KcT rj(j+i) 



■p(4) _ tv^Cp 



n^Pj^n'^Sij 2777, 

(5) Finally, the contribution of the operator — icf" [5^ cr] ■ (eegE*^™ x Vr)5'/(4m^) is 



(85) 






n3Pj-s-n'35i7 27772 



^(^+1) 2 



H,{nl ^ n'O) + 2Ii{nl ^ n'O) 



B. Quarkonium state corrections of relative order v"^ 

The quarkonium state (138 p is not an eigenstate of the complete Hamiltonian of pN- 
RQCD. The eigenstate may be constructed from fl55]) by systematically adding higher-order 
corrections, which are perturbative in the relative velocity v. Corrections may come from 
higher-order potentials {1/m and 1/m'^ terms) and from higher Fock states, which account 
for the coupling of the quark-antiquark singlet state to other low-energy degrees of freedom. 
We have to include such corrections at relative order v"^, both in the initial and in the final 
quarkonium states, in order to achieve a precision of relative order f ^ in the El transition 
rates. 

1. Corrections due to higher- order potentials 

The first-order correction to the quarkonium state fl38p induced by a correction 6hs to 
the Hamiltonian h^ is given by 

d^P' ^^ |/7'(P',A'))(°^ (°n^'(P'>A')| 



^ \ / UI-/-IJ 



771(0) _ p(0) 



X 



I £R f £r Tr {S^5/i5(r)S} \H{P, A))(°) . (87) 



We assume that, in order to account for corrections of relative order v"^, we need to include 
in Shs all the 1/m and l/m"^ potentials and, at order 1/m^, the first relativistic correction to 
the kinetic energy. Such a counting, which holds for weakly-coupled quarkonia, appears to 
be generally consistent with heavy quarkonium spectroscopy [16], and it is indeed the most 
widely used. It should be remarked, tawever. that m the case of strougly-eoupled quarkouia 
it is not the most conservative one [2S] . 
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In the assumed power counting, at relative order f^, Shs has the form 

Shsir) = -^ + SVsir) , (88) 

where p = —i'Vr- SVs is organized in powers of 1/m, 

where, at order 1/m^, we have distinguished between spin independent (SI) and spin depen- 
dent (SD) terms [ij, 

Vg^{r) =V}^\r) + l{V^?{r),p^} + ^L^ , (90) 

V^Sir) =vlS{r)L ■ S + V^l\r)S' + V^f^{r)S,,{r) , (91) 

with S = Si + S2 = (cTi + (T2)/2, L = r X p and Si2(f) = 3(f ■ cri)(r ■ (T2) — cr\ ■ (T2- Terms 
involving the center of mass momentum, which is — k^ for the quarkonium in the final state, 
are suppressed by an extra v and have been neglected. In the weak-coupling case, the above 



potentials read at leading (non- vanishing) order in perturbation theory (see e.g. 29j) 






fSM = ^=S?. t'i?w = 2;i^«^w, vA.'w = ^- m 



In the strong-coupling case, the potentials are non-perturbative and can be expressed in 



terms of Wilson loops to be eventually evaluated on the lattice 23j. 

There are, however, some observations that can be made without relying on any specific 
form of the potential, but just on its general structure ( l89l) -( l9T1) . Let us consider the radiative 
transition n^Pj — )• n'^Si 7. Initial state corrections due to 6hs read 



^n3P,^n'35.. = -^'Hn'^S.^ljd'RCE.ln'Pj) 



(1) 



_ ^ W(m3p,|^/.,(r)|n3p,)(°) („) 

2_^ (0) _ „(0) '^m3Pj^n'35i7' y^"^) 

m^n nl nil 

where E)^1 is the leading-order binding energy of a quarkonium with principal quantum 
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number n and orbital angular momentum L. Final state corrections read 

^ ^ ^<^Hn''Si\6hs{r)\m^Si)^<^^ (o) 

Z^ p(0) _ p(0) •^n3Pj^m3Si7 

"^ Z^ p(0) _ p{0) ■^n3Pj^m3Dj, 7 • W'*>' 

m,J' -'^n'O ^m2 

In both expressions, a sum over the intermediate state polarizations is understood. We have 
made use of the selection rule for the electric dipole matrix element and of the fact that the 
potentials flS^ -f pT]) do not change the total angular momentum of the state, which follows 
from [30] 

/r,25+lr |T/(2)/^\T Qlr>'25'+lr' \ _X X X J{J^^1)^L{L + 1) - S{S + 1) 

[n Lj\V^g[r) L ■ ^\n Lj,) — Oll'Ojj'Oss' 



2 

2S+lr I T/(2)/ \ \f2S+l 



x{n''-''Lj\V^'>{r)\n''''-^^Lj), (95) 
{n''^'Lj\V^f{r)S'\n''''^'L'j,) = 6ll'6jj'6ss' S{S + I) {n''^'Lj\V$\r)\n'''+'Lj), 

(96) 
{n''^'Lj\V^f^{r)SMr)W'^'Lj,) = 5jj,5ss'{S,2)Ljs{n''^'Lj\V^'lir) \n'''^'Lj), (97) 

where (812)101 = -4, (812)111 = 2, (812)121 = -2/5, (812)075 = and (8i2)ljo = 0. 
For the final state correction, we have assumed a completely non-degenerate spectrum. 
To complete the spin structure of the wave-function corrections, we need to compute 
^^^n^Si\ Vg_^!^{r) \m^Dj)^^\ It turns out that this matrix element does not vanish only for 
J = 1 and in this case 

POO 

^^\n^S^\ V^f^ir) |m=^Di)(°) = -2V2e^3n, ■ e^ss, / drr^ Rno{r)Rm2ir) VJ'^^(r) . (98) 

Jo 

Summing over the intermediate state polarizations allows to factorize the leading order 
amplitude. For instance, in the case of the contribution of the spin-tensor potential 
l/^J(r)8i2(f)/'m^ to the initial and final states, we obtain 



•^n3Pj-i.n'35i7 



/■oo 
m^n 111 ml ^ ' 

POO 

. / drr^RMr)Rn.2ir)V^^lir) 



-'n'O ^m2 
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Analogous results hold for the other potentials. The complete correction to the transition 
width coming from higher-order potentials is of the form 

■ph.o. potentials _ p5=l/ j\ (^nT\ 

J-„3pj^„/35j^ - -«-„„' WJ , l-LUiJ 

■ph.o. potentials _ p5=0 (^no^ 

where the general spin structure of R^^}{J) is R^^}{J) = Ann' + '^Bnn' + Cnn'J{J + 1) + 
-D„„'(Si2)iji {Bnn' is the correction coming from the spin-spin potential) and the general 
spin structure of Rnn''^ is -Rf^/° = Ann' + 2C„„'. The specific values of the coefficients Ann', 
Bnn', Cnn' and £*„„', which involve expressions similar to the ones displayed in Eqs. (1^ 
and fllOOp . depend on the form of the potentials and will not be discussed further in this 
work. 

2. Corrections due to higher Fock states 

The quarkonium initial and final states may also get corrections from the coupling of the 
heavy quark-antiquark pair to other low-energy degrees of freedom. We call these corrections 
higher Fock state corrections. 

For strongly-coupled quarkonia, we argued that we can neglect couplings with other low- 
energy degrees of freedom, see Sec. IIIBi In this case, we do not have new corrections coming 
from higher Fock states. 

For weakly-coupled quarkonia, we have to account for the coupling with low-energy glu- 
ons, see Eq. ([7]). A higher Fock state contributing at order v"^ is made of a gluon and a 
heavy quark-antiquark pair in a color-octet configuration. Moreover, photons may couple 
to a quark-antiquark octet state through fITOl) . Again we consider the radiative transition 
n^Pj — !■ n'^Si'j, whose relevant diagrams at relative order v"^ are shown in Fig. [HI 

The first two diagrams correspond to the normalization of the initial and final states, 
(SZ„3p, + 6Zn'3Sj^)/2. They can be calculated from the self energy of the state, first derived 
in 181], which is given for a generic quarkonium H by^ 

? r°° 



x(°)(/7(0, A)|r^e-'('^o'-<'V|iy(0, A))(°) , (103) 



With respect to [18|, a factor 1/Nc has been reabsorbed into the normahzation of the vacuum state. 
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FIG. 8: Color-octet contributions to the El transition n^Pj — )■ n'^Sij for weakly-coupled states. 
The double line stands for an intermediate octet state. 



where the Wilson line in the adjoint representation is 



^ -ig f dt'Ao{R,t')'^' 
e Jo 



V 



J 



(104) 



ab 



and hp = — V^/m + Vq. Deriving the self energy with respect to the energy provides the 



state normalization 



SZ 



dSE 



H{X) 



"^^^ - ^ I dtt{vac\gE'''{K,t)(j){t,0)'^^gE^\K,0)\vsic) 



OeP 6,0 



X ^^> {H{0, A) \r^e-'^''o'-^H)tr^H{0, X)y^> . (105) 



(0) 
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Thus, the contribution from the state normahzations to the amphtude reads 

X 






Diagram 2, the initial state correction, yields 

^^.^s., = -lY. :i^r^2V rfMvac|,E-(R,t)</,(t,0)-^,i?-(R,0)|vac) 

while diagrams 3a and 3b, the final state corrections, yield 

A%f::,.s.^ = -qY: :Z'Ztm / dt{.^c\gE-iK,mt,0):^W\K,0)\-^^) 

x<»>(n"'Si|r'e-l''"-<'i"r'|,n''Si)l°' , (108) 

^Fig.lll 3b vanishes, since r^ q~^^-^o ~-^,i'o'* r^ is a scalar that cannot change the angular momen- 
tum. The contribution from the last diagram, where the photon is coupled to an intermediate 
octet state, reads 

According to the power counting of Sec. IIIB[ all amplitudes f ll06p - (lll0p contribute to relative 
order l/{mv)^ x Aqqj^ ~ f ^ or l/{mv)'^ x Aq^q x l/[mv'^) ~ v^, where the factor l/{mv)'^ 
comes from the two r^ and we have assumed the chromoelectric correlator to scale like 

It is noteworthy that, in contrast to Ml transitions [l3|, weak-coupling non-perturbative 
contributions do not cancel for El transitions. This is due to the fact that the leading 
electric dipole operator does not commute with the kinetic energy. Hence, even in the weak 
coupling. El transitions are affected at order v"^ by non-perturbative contributions. These 
are encoded in the chromoelectric correlator {va,c\gE^'^(Il,t)(j){t,Q)^jj'^gE^^{'R,0)\vac). For 
E ^ Aqcd, the correlator reduces to the gluon condensate, which factorizes; however, this 
scale hierarchy is unlikely to be of relevance for El transitions. 
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VI. RESULTS AND COMPARISON WITH THE LITERATURE 



Now we can summarize our results and compare them with 10(]. The results are obtained 
by summing up all corrections calculated in the previous section, which include contributions 
coming from El operators of relative order f^, and initial and final state corrections. The 
complete decay rates n^Pj — )■ n''^Si'~f and n^Pj — )■ n'-^So'y, this last one obtained by leaving 
out spin-dependent contributions from the former expression, read to order k^/rn^ 



'J{J+1] 



<0) 



l + ^l'(^)-7±- 



k^ k"^^ I^{nl -^ n'O) 



Qm 60/3(^1 — )■ n'O) 



+ 



-'- n-i-Pi— +n"-5o 7 






h{nl^ n'O) 



.(0) 



1,05=0 ^7 ^Hnl^n'O) 



(111) 

(112) 



Qm 60 Iz{n\ — )■ n'O) 

Rnn'^iJ) and Rnn'^ are the initial and final state corrections; they include corrections coming 
from higher-order potentials, see Sec. IV B 11 and, in the case of weakly-coupled quarkonia, the 
color-octet contributions computed in Eqs. (ll06p -( fTT0l) . We have kept terms proportional to 
the anomalous magnetic moment, k^, for the sake of the following comparison, though these 
terms are suppressed by as{m) and go beyond our accuracy. Analogously, the expressions 
for the decay rates n^Si — )■ n'^Pj'y and n^So — ?■ n'^Pi'-f read 



2J + 1, 



n^^Sl—^n''^PJJ 



m 



l + fl-.(J) + ii *.Wl^nO) 



+ l^i^±i-2 



rU'-v 



(i + 4°)^ + :;3(i + 2«' 



Qm 60/3(n'l -)■ nO) 

1 .. ^ . emMn'l ^ nO) + 2h{n'l ^ nO) 

-(1 + 2/to ) 
m,- 



Isin'l ^nO) 



rnlSo-i-n'lPi7 



3r 



(0) 
nn' 



iS=0 



1 + Rr,„" + 



The decay rates derived in [lO| read 



k^ _ k^I^{n'l-^nO) 
Qm 60/3(n'l — )■ raO) 



(113) 
(114) 



lid 

n^Pj—^n'^Si 7 



n^Si— >-n'3Pj7 



.(0) 



1 + Rnn' (J) 



^ I^inl -^ n'O) k^ r^{nl -> n'O) k. 



+ 



^ 2m 



10 /3(nl ^ n'O) 2m I^inl -^ n'O) 
fc^ fJ{J+l) 



7 
m 



2J + 1. 



.(0) 



1 + -Rn.n' (•^) 



^ hin'l -^ nO) k^ r^{n'l -^ nO) 



(115) 



-Kr 



10 hin'l -^ nO) 2m Is{n'l -> nO) 
k^ fJ{J+l) 






2m 



-2 



(116) 



34 



The expressions appear different from tlie ones derived in tliis work because tlie basis of 
operators used in 10| is different from tlie one of pNRQCD. Tlie two basis are however 
related by a field redefinition so that at the end the two results are equivalent. To see this 
explicitly at the level of the transition widths, consider the radial Schrodinger equation 

\ ( d' 2 d\ L{L + l] 



EnlRnL 



d^ 2d 

m V dr'^ r dr 



rar^ 



+ Vl'\r) 



R. 



'TlL ■ 



From there, it follows (up to corrections of relative order v"^) 



k^h{nl -^ n 0) = — [6/;(nl -^ n 0) + 14/3(nl ^ n'O)] 

k^h{nl -^ n'O) = — [2/^(nl -> n'O) + 4Ji(nl ^ n'O)] , 

-k^h{n'l ^ nO) = — [6/^(n'l -^ nO) + 14/3(n'l ^ nO)] 

-k^hin'l ^ nO) = — [2/'(n'l -^ nO) + 4/i(n'l -> nO)] , 

m 



(117) 

(118) 
(119) 
(120) 
(121) 



where k^ = E^^ - ^Sl + C>{v^) in the first two lines and k^ = E^^^l - E^^,\ + 0{v^) in the 
third and fourth line. Using the identities (lll8p -f lT2Tl) . Eq. (Illip can be cast in the form of 
Eq. (IllSp and Eq. (I113p in the form of Eq. (I116p . Finally, a third way of presenting the 
same results, but in a more compact form, is 



n^Pj—m'^Si 7 



^(0) 



1 + Rnn' (J) 



^/5(nl-^n'0) k^ 
60 Isinl ^ n'O) ~ 6^ 



+^Q^ 



k^ fJ{J + l] 



em 7 

2m 



-2 



(122) 



2J + 1, 



n^Si—^n'^Pj 7 



^(0) 



1 + Rnn' (J) 



fc^J5(n'l^nO) k^ 
60 /3(n'l — )■ nO) 6m 

2m, 



— K 



h 

cm "^7 



-2 



(123) 



In summary, the decay width r„3pj_j.„/35j^ may be written up to order fc^/m^ in the 
equivalent ways (II lip . (IllSp or (I122p . the decay width Tn3Si^n'^Pj-y in the equivalent ways 
(I lisp . (I116P or (I123p . the decay width r„ip^_^„/i5g^ as in Eq. (11120 and the decay width 
rni5o->n'iPi7 as in Eq. (I114p . The obtained expressions are valid both for weakly-coupled 
and strongly-coupled quarkonia, the difference between the two cases being in the wave 
functions. Initial- and final-state wave functions affect crucially electric dipole transitions. 
The leading-order width, FJ^^,, depends on the wave functions, and at higher orders the 
integrals In and the initial- and final-state corrections R^^,^{J) and R^n'^ depend on the 
wave functions. 
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In the case of weakly-coupled quarkonia, the wave functions are Coulombic, which implies 
that r^^,, In and the initial- and final-state corrections due to higher-order potentials may 



be calculated in perturbation theory. The relevant potentials are those listed in Eqs. 
022]) • Weakly-coupled initial and final states get also corrections due to color-octet quark 
antiquark states coupled to low-energy gluons, which are parametrically of the same order as 
the other corrections. Color-octet corrections are given by Eqs. fll06l) - flll0p and depend on 
the correlator of two chromoelectric fields, which is a non-perturbative quantity. Therefore, 
at relative order v"^ even El transitions of weakly-coupled quarkonia are affected by non- 
perturbative corrections. 

In the case of strongly-coupled quarkonia, the potentials and, hence, the wave functions 
are non-Coulombic. The integrals Jjv and the wave-function corrections due to higher- 
order potentials, which are encoded in the coefficients Ann', Bnn', Cnn' and Dnn' defined in 
Sec. IV B 1\ are non-perturbative parameters. These non-perturbative parameters may be 
either derived from the quarkonium potentials evaluated on the lattice or fitted to the data. 
For strongly-coupled quarkonia, there are no relevant higher Fock state contributions to the 
initial and final states to be included. 

VII. CONCLUSIONS AND OUTLOOK 

The paper completes the analysis of radiative transitions in an EFT framework initiated 
in [13| with the study of Ml transitions. The EFTs are NRQCD and pNRQCD. 

The paper deals with El transitions, which are studied at relative order v"^, corresponding 
to order k^/m? in the transition width. All the relevant operators of pNRQCD are listed in 
Eq. (|9]). The matching, performed in Sec. IIIH shows that, if charm-loop effects are neglected, 
these operators do not get corrections from the momentum-transfer scale and keep the value 
inherited from NRQCD to all orders in perturbation theory and non-perturbatively. Charm- 
loop effects may be treated perturbatively and affect the matching beyond our accuracy. This 
non-obvious outcome may be considered the main result of the paper. 

As a consequence of the exact matching, we could provide the transition widths 
r„3p^^„/3sj^, r„35.j^„/3pj^, Tn^p^^n'^Soi ^"^^ ^ n^ So-^n' ^ p^ -f up to ordcr fc^/m^ for both 
weakly- and strongly-coupled quarkonia. Weakly-coupled quarkonia are those bound by a 
Coulombic potential, possible states being the quarkonium ground states; strongly-coupled 
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quarkonia are those bound by a non-perturbative potential, which eventually becomes con- 
fining in the long range. Strongly-coupled quarkonia are likely all states above the ground 
state. The transition widths have the same expressions for weakly- and strongly-coupled 
quarkonia, the only difi^erence lying in the wave functions and ultimately in the potentials. 
Weakly-coupled states also get corrections from intermediate quark-antiquark color octet 
states. The final expressions of the transition widths are listed in Sec. IVII Many alternative 
expressions for the widths are possible, all of them equivalent at order k^/rn^. We have listed 
some of them in the case of n^Pj — )■ n'^Si 7 and rc'Si -^ n'^Pj 7 transitions, equations (1122p 
and f ll23p providing the most compact expressions. 

The expressions for the widths agree, with some specifications, with the expressions 
obtained in [lOJ by reducing some covariant two-particle bound state equation. The specifi- 
cations are the foUowings. The expressions for the transition widths are valid up to relative 
order v"^. At this order the anomalous magnetic moment of the quark, k'q^, does not need 
to be included. If the anomalous magnetic moment is included, its expression is (j4]). This 
amounts to a small positive quantity of order as{m). No large non-perturbative correction 
affects Kg", as sometimes required in phenomenological treatments. According to a com- 
monly used power counting, wave-function corrections of relative order f ^ are induced by 
the potentials listed in Eqs. ( 189|) -( 19T|) . Typically, corrections induced by the 1/m poten- 
tial, V"*^^-*, have been neglected in the past, for the potential does not show up at tree level. 
For weakly-coupled quarkonia, non-perturbative corrections induced by low-energy gluons 
coupled to color octet quark-antiquark states have to be included as well. These correc- 
tions have been computed here for the first time and have not been included in any earlier 
treatment although they contribute at relative order v'^. 

Relativistic corrections to El transitions have, in some respects, opposite characteristics 
to the ones to Ml transitions. Allowed Ml transitions between quarkonium ground states can 
be described at relative order v"^ entirely in perturbation theory. In contrast. El transitions, 
even between weakly-coupled quarkonia, require at relative order v"^ a non-perturbative 
input. The reason is that, while the magnetic dipole operator commutes with the kinetic 
energy, leading eventually to the cancellation of the octet contributions, the electric dipole 
operator does not. Ml transitions between strongly-coupled quarkonia require at relative 
order v^ the non-perturbative matching of a yet unknown l/m^ operator. In contrast. El 
transitions between strongly-coupled quarkonia involve at most l/m,^ operators, which are 
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exactly known. Hence, a first principle calculation of El transitions at relative order v"^ 
is at present possible for all quarkonium states. Clearly, in the case of strongly coupled 
quarkonia, this requires parameterizing the lattice quarkonium potentials and solving the 
corresponding Schrodinger equation. 

Future applications of the present work include the numerical determination of the El 
transition widths between all S- and P-wave quarkonium states from the expressions given 
in Sec. IVII A consistent determination would require parameterizing the long distance be- 
haviour of the quarkonium potentials as evaluated on the lattice (for recent lattice results 



see, for instance. 
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ll-l33| ) . and matching it with the known short distance behaviour |l4l. |34|. 
solving the corresponding Schrodinger equation and finally evaluating the integrals I^ and 
the wave-functions corrections due to higher-order potentials. For weakly-coupled quarkonia, 
a parameterization of the chromoelectric field correlator would also be necessary. 

Finally, we mention that the EFT approach for quarkonium radiative transitions dis- 
cussed here can be translated to other systems beyond QCD. For instance, one could study 
atomic dipole transitions or dipole transitions in quirkonium, which is a candidate for dark 
matter 35|. The coupling constants of these systems are small, making them suited for a 
perturbative treatment. 
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Appendix A: Matching of Ml operators at 0{l/m^) 

Three 1/m^ operators contribute to Ml transitions at relative order v'^. One is the 
operator of Eq. f l5B]) . the other two are 

^ V^ {S\a-eeQB'^}VlS, (Al) 



and 

1 (Vr-o-)(Vr-S) , , . .- . . ^ ^ 

^-^V ^ {St,o-^eeQB^-^}V;V^,S. (A2) 

An even number of momenta, —i'Vr, is required by time-reversal symmetry; our analysis 
shows that other possible operators with two derivatives do not get contributions from the 
matching. 

♦ • * «N 



I I 







/ D- \ 



I I 

-i m. i •f ®- ^^^^-ft^ — 0^ — ®- 



FIG. 9: Some diagrams contributing to the matching at order 1/m^. 

We focus on the operators ( lAip and ( lA2p . They come from matching spin-dependent 
NRQCD amplitudes that contain operators with at least two spatial derivatives (see Fig. |9]). 
We have the following cases. 

1. Amplitudes that are made of a Oll/m^) operator of NRQCD may couple to longitu- 
dinal photons, to photons embedded in covariant derivatives or to the magnetic fields of the 
operators. The first case is excluded, for only transverse photons contribute to the operators 
(lAll) and (IA2I) . The second case does not contribute to the operators (1A1I) and (IA2I) . for 
spin-dependent 0{l/m^) operators of NRQCD contain at most two covariant derivatives. 
The third case involves the NRQCD operators 

6m Gin 

^^Hd2, a- ■ eegB'^jij - ^^J^Wa- ■ eeQW^^D'^ , (A3) 

^V^t (^ . D ^^^-Qera . Q + Q . ^^^-Qera ^ . D) ^^ , (A4) 
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and the corresponding c.c. ones, with c^^ — c^2 = 1 ^^^ c^™ = 0{as) [13|, l20|. For the 
purpose of matching the operators flAip and flA2p . the covariant derivatives may be replaced 
by simple derivatives that do not act on the gluon fields. We are, therefore, in the situation 
discussed in Sec. 1111 C[ paragraph (iv), and Eq. f l2U]) leads to 



Vr "B (Vr-cr) (Vr'-B) 



V^^ = 1, V ^ = cT . (A5) 



These relations are exact because, as we are going to detail in the following, NRQCD am- 
plitudes with insertions of 1/m and 1/m^ operators factorize at relative order v'^ into con- 
tributions to 1/m and 1/m^ pNRQCD operators. 

2. In diagrams made of a kinetic energy and a chromoelectric spin-orbit operator, the 
photon is coupled to one of the covariant derivatives. If it is coupled to the covariant 
derivative of the chromoelectric spin-orbit operator, then the two derivatives come from the 
kinetic energy operator, which factorizes. If it is coupled to the kinetic energy operator, 
then the photon field multiplies a derivative, the operator commutes with the gluons and 
factorizes. 

3. The situation is similar for diagrams with a kinetic energy and an electric spin-orbit 
operator insertion. 

4. In diagrams made of two kinetic energy and a chromomagnetic dipole operator, the 
photon is coupled to one of the kinetic energy operators. If the photon field multiplies a 
derivative, then either the operator commutes with the gluons and factorizes or it does not 
and then the two derivatives come from the other kinetic energy operator, which factorizes. 
If the photon field multiplies a gluon field then the two derivatives come from the other 
kinetic energy operator, which factorizes. 

5. The kinetic energy operator also factorizes in diagrams with a magnetic and a chro- 
momagnetic dipole operator. 

6. Finally, in diagrams made of two kinetic energy and a magnetic dipole operator, 
the magnetic dipole operator commutes at relative order v"^ with the kinetic energy and 
factorizes. 

We conclude therefore that ( JA5I) holds to any order of the matching, while only the Wilson 

ctS 

coefficient V^ of ( 1361) . the only operator without derivatives, may possibly get corrections 
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from the matching beyond tree leveL 
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